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Abstract Gravitationally lensed QSOs (GLQs) at 1 
< z < 2 play a key role in understanding the cos- 
mic evolution of the innermost parts of active galaxies 
(black holes, accretion disks, coronas and internal jets), 
as well as the structure of galaxies at intermediate red- 
shifts. With respect to studies of normal QSOs, GLQ 
programmes have several advantages. For example, a 
monitoring of GLQs may lead to unambiguous detec- 
tions of intrinsic and extrinsic variations. Both kinds of 
variations can be used to discuss central engines in dis- 
tant QSOs, and mass distributions and compositions of 
lensing galaxies. In this context, UV data are of partic- 
ular interest, since they correspond to emissions from 
the immediate surroundings of the supermassive black 
hole. We describe some observation strategies to anal- 
yse optically bright GLQs at z ~ 1.5, using ISSIS (CfS) 
on board World Space Observatory-Ultraviolet. 

Keywords Gravitational lensing; Quasars: black holes 
and accretion; Galaxies: structure 



1 Gravitationally lensed QSOs in the UV: 
central engines 

Clusters, groups and individual galaxies are so massive 
and compact that light rays from background sources 
are gravitationally deflected by them. In the strong 

L. J. Goicoechea 

Departamento de Ffsica Moderna, Universidad de Cantabria, 
Avda. de Los Castros s/n, 39005 Santander, Spain 

V. N. Shalyapin 

Institute for Radiophysics and Electronics, National Academy of 
Sciences of Ukraine, 12 Proskura St., 61085 Kharkov, Ukraine 

R. Gil-Merino 

Departamento de Ffsica Moderna, Universidad de Cantabria, 
Avda. de Los Castros s/n, 39005 Santander, Spain 



lensing regime, the involved gravitational fields bend 
light to form several images of the same background 
objec t ( Blandford fc Naravan 19921 Schneider et al 



19921: IWambseanss lll998HSchneider etal. II2006I ). The 

fir st gravitationally len sed QSO (GLQ) was discovered 
bv lWalsh et al. 1 dl979h 30 years ago, and there are cur- 
rently more than 100 known GLQs. Most of them are 



active galactic nuclei (AGN ) at z > 1 (jKochanek et al. 
20ldllnada fc Oeuri ilioioh . 



We initially consider the population of optically 
bright GLQs. From the CASTLES^ and SQLS3 
databases, we select 46 GLQs showing at least two 
images with V < 20 mag. This GLQ sample is de- 
picted in Fig. [T] (only two sample members with angu- 
lar size A9 > 10" are not included in the plot, where 
A9 represents the maximum separation of any pair of 
images, twice the average distance of images from the 
lens center or something like that). We also classify the 
population in terms of the U-band flux of the brightest 
image of each GLQ (V*): V* < 17.5 mag (blue circles), 
17.5 < "V* < 18.5 mag (green squares), and V* > 18.5 
mag (red triangles). About 50% of the GLQs with V* < 
17.5 mag have their redshifts in the interval 1 < z < 
2 (see the region between both vertical dashed lines in 
Fig- QJ ! and we focus on this densely populated cosmic 
shell containing 61% (28 out of 46) of the whole CAS- 
TLES + SQLS sample, i.e., the population of bright 
GLQs at z ~ 1.5. 

A lensing galaxy (or group/cluster) acts as a natural 
telescope, so it produces a GLQ with magnified images 
(with respect to the single image of a non-lensed QSO 
at the same z and having the same intrinsic luminos- 
ity). However, the dust in the galaxy may significantly 
demagnify some GLQ images at certain wavelengths 

1 http://www.cfa.harvard.edu/castles/ 
; http://www-utap.p hys.s.u-tokyo.ac.jp/~sdss/sqls/| 



2 



CD 

O 
W 

u 



CD 

< 



CV2 



1 1 

— 

— 


1 1 1 1 
• 




i i 

- 


■ 




• A 






* 












• 


1 1 


i i i i 


i i i i 1 i 





4 



Fig. 1 CASTLES + SQLS sample of bright GLQs. We 
display angular sizes and redshifts for GLQs having at least 
two images with < 20 mag. Taking the V-band flux 
of the brightest image of each GLQ (V*) into account, the 
sample is divided in three different subsamples: V* < 17.5 
mag (blue circles), 17.5 < V* < 18.5 mag (green squares), 
and V* > 18.5 mag (red triangles). Two vertical lines define 
the region of interest (1 < z < 2) 



Falco et al. I (|1999l ). Moreover, the presence of high 
column densities of neutral hydrogen (in intervening 
objects at z < 1.5) could greatly hinder or make im- 
possible the detection of photons at \ b s < 2280 A, i.e., 
A < 912 A (= Lyman limit) in the GLQ rest frame 
[z ~ 1.5). At the redshift of interest, we remark that 
optical observations in the BVR Bessell and gr SDSS 
bands correspond to far/middle UV (FUV/MUV) emis- 
sion (A - 1700-2600 A). Hence, even in the worst case 
(total absorption below the Lyman limit), UV observa- 
tions at A o6s - 2300-4000 A (EUV/FUV emission: A ~ 
920-1600 A) are possible. These have an extraordinary 
relevance in the study of the innermost accretion flow 
in a distant QSO, and thus, to reveal its hottest and 
most compac t sources aroun d the supermassive black 



hole (SMBH: [Peterson Il2009h . 



2 Time domain studies of optically bright 
GLQs atz~ 1.5 

2.1 Intrinsic variability 

A monitoring campaign of a normal (non-lensed) QSO 
at z ~ 1.5 produces light curves in the UV/optical 
continuum whose variations are difficult to interpret: 
intrinsic variability or fluctuations caused by the in- 
tervening medium? However, the presence of a grav- 
itational deflector at intermediate redshift may result 
in several (usually two or four) images of the same 
QSO. As the ray paths are different for different im- 
ages, the corresponding magnifications and travcltimcs 



do not agree with each other. Thus, in the absence 
of significant extrinsic variations (see the opposite case 
below), each intrinsic flare is observed in two or more 
images with different arrival times. This allows ob- 
servers to fairly interpret the nature of the variability, 
and even more, the light curves of the images can be 
com bined to make one better-sampled intrinsic recor d 
(e.g. Kawaguchi et al. Ill998t Goicoechea et al. 1l2008a ). 

A structure function analysis is a powerful tech- 
nique to quantify typical luminosity variabilities at 
different rest frame lags and to check for tim e scales 



of intrinsic flares 



(e.g. 
l200lh 



Cid Fernandes et al. 



20001 : 



Collier fc Peterson I l200Lh . For two optically bright 
GLQs at z = 1.4 (SBS 0909+532 and Q0957+561), 
we recently reported the possible existence of time- 
symmetric flares (with symmetric r ise and decay) at 
A ~ 19 00-2600 A and lasting - 100 d (|Goicoechea et al. | 
2010al) . As both GLQs have similar SMBH mass 



( Peng et al. 1120061 ) . the symmetry and timescale could 
be related to this physical quantity. There is also evi- 
dence that 170-d asymmetric flares are generated in the 
nuclear region of Q0957+561, but these events seem 



to hav e an intermittent character (jGoicoechea et al. 
20081)). Future studies involving larger samples and 
longer observation periods will confirm/reject the 
possible correlation between symmetry/timescale and 
SMBH mass, as well as the generation of asymmetric 
flares in an intermittent way in some GLQs. The radio- 
loud GLQ Q0957+561 also has an achromatic structure 
function at A ~ 1900-2600 A, suggesting the existence 
of an illuminated accretion disk (reproces sing scenario; 
ICollin-Souffrin Ill99ll Ikrolik et al. Ill99ll ): high-energy 
(X-ray/EUV) flares from a source above the SMBH are 
reprocessed by disk rings into lower energy events. 

For a given AGN, time delays between variations at 
different UV/optical continuum wavelengths constrain 
the SMBH accretion physics. The standard reprocess- 
ing scenario consists of a very compact X-ray source 
located at a height of a few Schwarzschild radii (h ~ 
a few Rschw) above the SMBH, which illuminates a 
standard g eometrically-thin and optic ally-thick accre- 



tion disk (jShakura fc Sunvaev 1 119731 ) . This corona 



disk coupling scenario predicts light travel time de- 
lays between an X-ray central flare and its replicas (re- 
processed versions) in different rings of the accretion 
disk. The smallest radii correspond to the hottest re- 
gions and shortest wavelength emissions, so light travel 
time delays r = R/c (R > h is an arbitrary ra- 
dius and c is the speed of light) can b e rewritten as 



wave l ength-dependent lags r oc A 4//3 (e.g.lCollier et al 



19991 : ISergeev et al. 1120051: ICackett et al. Il2007t ). Light 



curves of many local AGN show UV/optical inter- 
wavelength lags supporting the standard reprocessing 



3 



scenario (e.g. Sergeev et al. I 2005 ). However, there are 
also indications that UV/optical fluctuations in several 
local AGN cannot be produced from standard repro- 
cessing of X-ray flares associated with observed X-ray 
variations. 
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Feb Mar Apr May Jun 



For example. I Arevalo et al. I ([20081 ) simulated thermally- 
reprocessed .B-band light curve s of the local QSO MR 
2251-1 78, using the technique by Kazanas fc Navakshin ^ 
( 20011 ) and the observed X-ray record as input. They 
showed that even observed short-timescale (lasting tens 
of days) B-band fluctuations do not agree with standard 
reprocessing simulations. If we concentrate on the two 
consecutive fluctuation s in the 3950-4100 (M JD-50000) 
period (see Fig. 9 in I Arevalo et al. 1 120081 ). only the 
green line seems to account for the B-band data. This 
line comes from a non-standard reprocessing simulation 
with h — 50 Rschw, i.e., the irradiating source does not 
correspond to a X-ray corona just above the SMBH. Al- 
though the X-ray source might be really located at the 
base of a jet, we note that at least one best-fit parame- 
ter, the observer latitude, is physically unplausible (an 
edge-on disk generates zero flux a t the observer and it is 
inconsistent with a type 1 AGN) . lArevalo et al. ( 20091 ) 
also showed that the reprocessor of rapid X-ray flares 
in the Seyfert nucleus NGC 3783 could correspond to 
the broad emission-line region (instead the accretion 
disk). In brief, short-timescale UV/optical variations 
are very likely related to reprocessing of high-energy 
events, but there is no standard scenario working in 
all local AGN. The main high-energy source may not 
be the X-ray corona (X rays from the base of a jet, 
EUV emission from a non-flat innermost region of the 
accretion flow, etc) and/or the reprocessor may be dif- 
ferent than a standard accretion disk (non-standard 
disk, broad emission- line region, dust torus, etc). 

Time lags between variations at different wave- 
lengths has been well studied only for local AGN. 
However, very recently we have detected a promi- 
nent optical variation in the image A of the GLQ 
Q0957+561, and searched for the imprints of a hy- 
pothetical intrinsic flare in the light curves of the B 
image, but 14 months lat er (the well-known time delay 
between both images, e.g. Shalvapin et al. Il2008l) . For- 
tunately, we have observed the GLQ (A6 ~ 6") with 
several facilities at several wavelengths: X-rays (Chan- 
dra), NUV (Swift/UVOT) and optical/NIR (Liverpool 
Robotic Telescope). In Fig. [2j we show preliminary 
NUV/optical/NIR light curves of Q0957+561B (until 
mid-May 2010) from our monitoring campaign with 
the UVOT on board Swift (black circles) and the Liv- 
erpool Robotic Telescope (rest of filled symbols) . Open 
symbols trace the time evolution of Q0957+561A (14 
months before) and demonstrate the intrinsic origin of 
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Fig. 2 Preliminary light curves of Q0957+561B over the 
first semester of 2010. Records of Q0957+561A (but 14 
months before) are also depicted for comparison purposes 
(open symbols; see main text). The UV/optical/NIR vari- 
ability of the B image of the distant GLQ has been monitor- 
ized with Swift/UVOT (black circles; the rest frame wave- 
length is A ~ 1440 A), and the Liverpool Robotic Telescope: 
A - 1940 A (filled blue squares), A ~ 2560 A (filled geen 
triangles), A ~ 3100 A (purple inverted triangles), and A ~ 
3710 A (red stars). The light curves are shifted in magni- 
tude to overlap them at day 5250 



the variability (e.g. the filled blue squares correspond 
to Q0957+561B data in the g SDSS band, and these 
are consistent with the open blue squares describing the 
behaviour of Q0957+561A in the same optical filter). 
Our NUV/optical/NIR follow-up covers the rest- frame 
spectral region A ~ 1440-3700 A, so we have lost EUV 
emissions at A ~ 920-1200 A, which probably come 
from the innermost accretion flow. In any case, the 
simple use of the U filter in the UVOT filter wheel al- 
lows us to properly analyse the FUV emission at A ~ 
1440 (8.6 eV). This is also expected to be created in 
internal rings of the accretion disk around the SMBH. 
At present, we are comparing the flux fluctuations of 
Q0957+561B from X-rays to NIR waves. 

Time-delayed responses of emission lines to contin- 
uum variations can be used to infer sizes of broad 
emission-line regions, while FWHMs of emission lines 
are directly related to the kinematics of line-emitting 
clouds. Both the geometry and kinematics of broad 
emission-line regions a re useful tools to estimate masses 
of central SMBHs fe.g. |Peterson~ll2009l ). Unfortunately, 
there are no spectral monitoring campaigns of GLQs 
(with the exception of COSMOGRAIL's recent efforts 
to follow-up the spectral variability of Q2237+0305 - 



•1 



1-3 



< 

LO 
O 

CO 

o 

+ 

co 

CV2 
CV2 



LO 
O 



~i 1 1 1 1 1 1 1 r 




4 



J I I I I I L_ 



_L 



1500 



1550 



Fig. 3 Einstein Cross at optical wavelengths. Apart from 
the four images A-D of the distant GLQ Q2237+0305, one 
can also see the local face-on spiral galaxy G2237+0305 act- 
ing as a gravitational lens 



see below, and an ongoing programme to take spectra 
of the two images of Q0957+561 on a regular basis), 
which might be much more efficients than programmes 
with non-lensed QSOs. This is because effective sam- 
pling periods are twice (double QSO) or even four times 
(quadruple QSO). 

Apart from unambiguous detections of intrinsic fluc- 
tuations and important improvements in effective sam- 
plings, it is clear that studies of GLQs have much more 
potential than those from non-lensed QSOs. Thus, 
time delays between lensed images are related to lens- 
ing mass distributions, and provide extremely valuable 
information on dark mat ter halos around galaxie s at in- 
termediate redshifts (e.g. Schneider et al. Il2006l ). From 



very distant GLQs, one can also infer unknown redshifts 
(not only dark halos) of lensing objects. For exam- 
ple, the time delays in the quadruple QSO H1413+117 
(Cloverleaf) at z = 2.56 have been used to refine the 
lens mod el and to estimate the previously unknown lens 
redshift ( Goicoechea k, Shalvapin IboiOblb 



2.2 Extrinsic variability 

Apart from the constant macrolens magnification (due 
to the lensing structure as a whole) of an image of a 
GLQ, an additional variable microlens magnification 
may be generated by compact objects (moving stars 
or MACHOs) in the main lensing galaxy (gravitational 



MJD-50000 

Fig. 4 GLITP V-band fluxes (mjy) of Q2237+0305A 
(blue circles) and Q2237+0305D (green triangles) between 
October 1999 and February 2000. The observations were 
made with the 2.56-m Nordic Optical Telescope, and the 
two light curves are shifted in flux to overlap them at day 
1450 



microlensing; e.g. Schneider et al. 20061) . In princi- 
ple, moving clouds of dust may also produce a vari- 



able transmission (e.g. iGoicoechea et al. 1 120051 ) . and 



thus, some extrinsic variability. However, most stud- 
ies focused on the microlensing interpretation of ex- 
trinsic signals in GLQs. When the extrinsic variability 
can be unambiguously sep arated from the intrinsic one 



(e.g. IParaficz et al~ll2006t . resulting extrinsic fluctua- 



tions may be related to the nature of the lensed and mi- 
crolensed source, as well as the composition of the main 



lensing galaxy (e.g. IShalvapin et al. Il2002[ iKochanek 



2004; Paraficz et al. 



200 



6) 



Microlensing (extrinsic) variability was initially dis- 
covered in the G LQ Q2237+0305 (Einstein Cross; 
Irwin et al. 1989). and we mainly comment results on 
this emblematic QSO. The gravitational lens system is 
formed by a distant QSO at z = 1.695 and a local lens 
galaxy (face-on spiral at z =0.039) that q uadruples the 



images of the QSO (jHuchra et al. II1985I ). This com 



pact lens system with A9 ~ 2" appears in Fig. [3] Some 
optical continuum monitoring campaigns have led to 
very accurate and reasonably well-sampled light curves 



(Wozniak et al. 


2000h. GLITP (Alcalde et al. 


2002) 


and COSMOGRAIL ( 


Eieenbrod et al. 20081) records 



The OGLE y-band records cover a time frame of about 
10 years, while the GLITP Fi?-band curves only last a 
few months but with a daily sampling, bad weather or 
technical problems as ide (see also the Ma idanak light 



curves in the R band; IVakulik et al. 1120061 ) . The COS 



MOGRAIL monitoring is the only spectroscopic follow- 
up (A > 1450 A in the rest frame of the QSO) of a GLQ 
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to date. In Fig. gj the GLITP U-band fluxes (mJy) of 
Q2237+0305A and Q2237+0305D arc plotted together 
as blue circles and green triangles, respectively. The 
time delay between the im ages A and D is less than 2 



days ( Vakulik et al. 1 120061 ) . so the time-axis in Fig. @] 



practically traces the emission time. There is a signifi- 
cant fluctuation in A, which is not detected in D at simi- 
lar emission times. This extrinsic fluctuation allowed us 
to discuss the U-band source intensity profile and size, 
the mass of the SMBH, the mass accretion rate, the 
lens transverse velocity, and the typi cal mass of stars 



in the bulge of the local spiral galaxy (IShalvapin et al. 



2002HGoicoechea et al. 1120031: iGil-Merino et al. II2006I) . 



The timescale of microlcnsing variability depends on 
the redshift of (distance to) the main lensing galaxy: 
more distant lenses produce longer variations. There- 
fore, the Einstein Cross is of particular interest because 
of the observed short-timescale microlensing fluctua- 
tions (e.g. Fig. [4]). For non-local lenses, one observes 
long-term microlensing effects that induce small (op- 
tical ) flux__variatipjiB^ 



(e.g. [ Gaynullina et al. 11200a iFohlmeister et al. 112007 



IShalvapin et al. 20091) . QSO emission regions have 



different sizes in different wavelengths, and microlens- 
ing effects are stronger f o r mor e compact sources. For 



example, iPoolev et al. I (|2007l) used observed X-ray 
and optical flux ratios between images of ten quadru- 
ple GLQs to show that X-ray sources (hot coronas 
above SMBHs?; see above) are much more strongly mi- 
crolensed (more compact) than sources associated with 
optical observations (accretion disks around SMBHs are 
prime candidates). Thus, the UV extrinsic variability 
should be stronger than the extrinsic variability at opti- 
cal wavelengths, and it should probe innermost regions 
of accretion disks. 



3 GLQ programmes using ISSIS (CfS) on 
board World Space Observatory-Ultraviolet 
(WSO-UV) 

The aim of the WSO-UV missio n is to study the Un i- 
verse in the UV spectral range ( Shustov et al. I [2009 ). 
but also including complementary optical/NIR capabil- 
ities (4000-10000 A). This I.70-m space telescope will 
be launched in 2014, and it will be the replacement 
for the Hubble Space Telescope (HST) in many ways. 
While the 4000-10000 A interval is not a problem for 
ground-based instruments, the 1000-4000 A range is 
problematic or beyond the sensitivity of ground-based 
telescopes. Therefore, WSO-UV has its main poten- 
tial in the UV wavelength interval. The mission is an 



Table 1 Exposure times required to obtain SNR = 100 (V 
— 18 mag) 



Filter 


Emission 


Time (s) 


F170W 


EUV 


- 430 


F255W 


EUV 


- 190 


F336W 


FUV 


- 25 



international collaboration, with Russi^, Spair0 and 
GermanjH playing the leading role. 

WSO-UV is equipped with multipurpose instrumen- 
tation to carry out spectroscopy and imaging. In partic- 
ular, we concentrate on the Channel for Surveys (CfS) 
of the Imaging and Slitless Spectrograph Instrument for 
Surveys (ISSIS). The ISSIS CCD (CfS) is characterized 
by 

1. small pixels (0.05"), 

2. narrow distributions of brightness for point-like 
sources (FWHM = 0.1"), and 

3. a field-of-view (FOV) = 3.4' x 3.4'. 

Hence, the CfS is ideally suited for GLQ programmes, 
which require a high angular resolution (see Fig. [T]) 
and a reasonably large FOV (to find reference stars). 
Its UV spectral coverage (1100-4000 A) corresponds to 
A ~ 440-1600 A in the rest-frame of bright GLQs at 
z ~ 1.5, and thus, the CfS can be used to probe the 
EUV-FUV emission of GLQs. We propose to follow-up 
the variability of several targets over observation peri- 
ods of about six months and using a dense sampling 
(e.g. two or three times per week). To optimise the 
monitoring of each GLQ, starting times (of observa- 
tion periods) would be decided through an optical alert 
system (just after detections of prominent optical gra- 
dients). As we remark in Sect. 2, this time-domain 
programme has high potential for different cosmolog- 
ical studies. Moreover, very preliminary simulations 
(using the ISSIS Exposure Time Calculator, Alpha 2 
version) suggest that it is not too time consuming. In 
Table [I] we present some exposure times (broad-band 
filters) that are needed to obtain a signal-to-noise ratio 
(SNR) of 100 for a target with V = 18 mag. Taking 
these (preliminary!) results into account, a whole cam- 
paign for a typical GLQ might last ~ 645 s/d x 50 d 
~ 30 ks (science time). 

We also propose a second programme: EUV database 
of GLQs. Its main scientific goals are: 

" http://wso.inasan.ru/ 

4 http: / /www. wso-uv.es/index.php/home. html 
^http://www.uni- tuebingen,de/index.php?i d=8711| 



6 



Table 2 Exposure times required to obtain SNR = 100 (V 
= 17-20 mag) 



Filter 


Time/F=17 (ks) 


Timc/U=20 (ks) 


F170W 


- 0.17 


- 2.63 


F255W 


- 0.08 


- 1.19 


F280N 


- 0.40 


~ 6.14 



To study nuclear and circumnuclear EUV emissions 
(A < 1200 A within a 0.3" radiu s) of a sleeted sample 
of GLQs Ce.g. lHutchings1l2003l) . 



• To compare with other observations (e.g. VLBI radio 
jets), and to contribute to a database covering the 
whole electromagnetic spectrum. 

• To reconstruct the EUV morphology of central re- 
gions in active galaxies, and to constraint gravita- 
tional lens scenarios. 

Table [2] shows exposure times required to obtain SNR 
= 100 for targets with 17 < V < 20 mag. The faintest 
targets (V = 20 mag) may only be reasonably well im- 
aged in a narrow filter from relatively long exposures 
(exceeding 6 ks). See the GLENDAMA web site at 



http://grupos.unican.es/glendama/index.htm for up- 



dates on our GLQ project with WSO-UV. 
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